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ABSTRACT 

From considerations of spin evolution and kinematics in the Galactic potential, we argue that the pulsars 
B1913+16, B1534+12, and B2127+11C may be younger than previously assumed, and we find that a lower 
bound on the formation and merger rate of close double neutron star binaries is lO -6 ' 7 /^ 1 /^ 1 y r ^ where fi, is the 
beaming fraction and /</ accounts for the possibility that the known NS-NS binaries are atypical of the underlying 
population (e.g., if most such binaries are born with shorter orbital periods). If we assume no prior knowledge of 
the detectable lifetimes of such systems, the rate could be as large as ~ 10~ 50 //~ '/7 1 Y 1 " 1 - From both plausible 
bounds on ft, and fj, and a revision of the independently derived limit proposed by Bailes (1996), we argue that 
a firm upper bound is 10~ 4 yr" 1 . We also present a unifying empirical overview of the spin-up of massive binary 
pulsars based on their distribution in spin period P and spin-down rate P, finding evidence for two distinct spin-up 
processes, one dominated by disk accretion, the other by wind accretion. We argue that the positions of binary 
pulsars in the P-P diagram can be understood if (1) there exists a Galactic population of pulsars in double neutron 
star systems with combinations of spin and orbital periods that have prevented their detection in surveys to date; 
(2) recycled pulsars in wide-orbit binaries are not born near the canonical spin-up line in the P-P diagram because 
they were predominantly spun up through wind accretion; and (3) there exists a disfavored evolutionary endpoint 
for radio pulsars with spin periods 5-30 ms and P > 10~ 19 s-s" 1 . 

Subject headings: binaries: close — gravitation — pulsars: individual (B1534+12, B1913+16) 



1. INTRODUCTION 

In anticipation of the new observational window afforded 
by interferometric gravitational-wave observatories now under 
construction, there has been considerable interest in the rates 
of formation and merger of binary systems containing at least 
one neutron star. For double neutron-star (NS-NS) systems, es- 
timates of these rates depend on the assumed "detectable life- 
times" of the radio pulsars in known close NS-NS binaries. 

The detectable lifetime of a known NS-NS binary is, presum- 
ably, the sum of the post-accretion age of the observed radio 
pulsar and the remaining timespan during which it will con- 
tinue to be observable. For the latter, some past authors have 
used the time to coalescence of the system, while others have 
argued that luminosity evolution limits the detectable lifetime. 
For age estimates, past work has relied exclusively on the "char- 
acteristic" spin-down timescale r c = P/2P, but r e is an accurate 
chronological age only if the present-day spin period P is much 
greater than the pulsar's period Pq at "birth" (i.e., at the end of 
the accretion phase for recycled pulsars), and for dipolar radia- 
tion braking with a constant magnetic field. Since Po may not 
be small compared to P, the spindown timescale is generally an 
overestimate of the time since spinup; consequently, birth and 
merger rates for NS-NS systems based on r e may be underesti- 
mates. 

We argue that a neutron star's rotational history since the ces- 
sation of spin-up may be used to better estimate its true age. 
For each known system, we check such age estimates against 
kinematic constraints derived from motion in the galactic grav- 
itational potential (Sections 2 and 3). In Section 0, we use our 
results to estimate the detectable lifetimes, and hence birthrates, 
of the known NS-NS binaries, revisiting the assumptions un- 
derlying earlier work. In addition, we show that there are kine- 
matic solutions for both B1534+12 and B1913+16 for which 



they are young (ages ~ 10 Myr); we describe how these ages, 
if characteristic of binary neutron stars, might alter estimates 
of their birth and merger rates. Finally, we show in Section ^ 
that the phenomena of spiral-in and merger of NS-NS systems 
are reflected in the P-P diagram, and argue that details of the 
spin-up process for other binary pulsars similarly leave behind 
important imprints. 

2. THE AGES OF NS-NS BINARIES 

2.1. The Spin-up Lines 

Through accretion of mass and angular momentum from 
an evolving companion star, a neutron star can be spun up 
to an equilibrium period approximated by the Kepler orbital 
period at the Alfven radius (Ghosh & Lamb 1992), P eq = 

0.95 sec £ 3 72 /^3q 7 m" 3 / 7 kJIr^ 1 M l ^ , where fi^o is the neu- 
tron star's magnetic moment expressed in units of 10 30 Gem 3 , 
fto.4 is the opacity of the accreting material in units of 0.4 
cm 2 g -1 (the Thomson opacity for ionized hydrogen), is the 
stellar radius in units of 10 6 cm, Mi .4 is the neutron star's mass 
in units of 1.4 M Q , and m is the mass accretion rate in units of 
the Eddington rate (accretion onto a sphere is assumed). The 
quantity £ = ^/ui 2 J^ is of order unity, where f ~ 0.5 encapsu- 
lates uncertainty in the Alfven radius and oj s < 1 is the fastness 
parameter. In the magnetic dipole model, PP = 87r 2 ?7/i 2 /(3/c 3 ), 
where / is the stellar moment of inertia and 77 <J 1 is a geomet- 
rical factor, which then implies the "spin-up line" relationship 
between P and P for the endpoint of spin-up by accretion, 

P=aP 4/i , (1) 
a = 1 . 10 x 10 _ 15 s" 4/3 77 C" 7/2 m n~ Q \ R 6 m\{1 1^ , (2) 

where P is in seconds, and / = Z15 10 45 g cm 2 . 
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FIG. 1. — A portion of the P-P diagram for pulsars in the Galactic disk 
and B2127+11C in the M15 cluster. Pulsars are shown as dots, with 
encircled dots denoting binary systems; the two empty circles repre- 
sent the eclipsing binaries. The sizes of the circles scale as logP or (,; 
heavy circles denote companion masses greater than 0.45 Mq, i.e., 
likely CO rather than He white-dwarf companions, and ellipses repre- 
sent double neutron-star binaries. Spin-down rates have been corrected 
for apparent acceleration due to proper motion and differential galactic 
acceleration (e.g., Camilo, Thorsett & Kulkarni 1994) where such cor- 
rections are known. For PSRs B 1534+12, B1913+16, and B2127+1 1C 
(M15C) evolutionary tracks for rotational spin-down assuming dipolar 
magnetic braking are shown — tick marks delimit time intervals of 20 
Myr; smaller ticks indicate the portion of each track for which P m ±, < 4 
hours. The tracks end at the expected coalescence time for each sys- 
tem, having begun on the equilibrium spin-up line (eq. |l|), shown as a 
solid line within a gray region; the latter indicates the extent to which 
the position of this line may vary if assumed quantities (such as mass 
accretion rate) take on somewhat different values. The spin-up line for 
face < feq, equation (^), assumes f acc = 10 6 yr. The Hubble limit is the 
locus of points for which P/(2P) = 20 Gyr, and the death line is the 
pair-creation limit signifying turn-off of radio emission. 

While the scaling P oc P 4 / 3 in equation ([[]) is well-founded 
for neutron stars with predominantly dipolar magnetic fields, 
there are numerous sources of uncertainty in the coefficient 
a: the parameters £, ui s , and r\ are not well known, and sub- 
Eddington accretion, m < 1, is likely (e.g., if inflow onto the 
neutron star is restricted to the polar caps). Also, the preferred 
spin-up scenario for NS-NS systems, accretion from the helium 
core of a companion stripped of its hydrogen envelope, suggests 
an opacity kqa = 0.5. We adopt a "fiducial" spin-up relation 
with coefficient a= 1.1 x 10~ 15 s" 4 / 3 and perform calculations 
for a range of values. 

If the evolutionary timescale of the donor star is short, the 
accretion phase may not last long enough to spin the neu- 
tron star up to equilibrium. The time to reach equilibrium 
is f eq ~ I/(Mr\); if mass accretion lasts a time f acc = 10 6 f acCl 6 
yr < the spin period at the end of the accretion phase is 
P = ^eqfeq Aacc > Peq, which implies the relationship 

P P? = 3.60 x 10- 20 ry C 1 ' 2 cuf m 6 C 6 4 4a R~e M? A , (3) 

where we have used P = lOPio ms. Equation (|j) is extremely 
sensitive to all of the parameters, and represents a family of 



P oc P spin-up lines in the P-P diagram; the line for m - 1 and 
?eq = 10 6 yr is plotted in Fig. 1. The intersection of this alternate 
spin-up line with the equilibrium P 4 / 3 line defines a minimum 
post-accretion period (corresponding to f acc = f eq ) that depends 
on the total mass accreted, but is insensitive to other parameters 
(see also Kulkarni & Phinney 1994). If we include the effects of 
accretion-induced field decay following Shibazaki et al. (1989), 
we find no qualitative change in the form of the non-equilibrium 
spin-up line. 

2.2. Spin-down Ages 

The time to spin down from birth period Pq to present-day 
period P, with constant braking index n, is 



T=T C [i-(p /pr 1 ] 



(4) 



where the characteristic age t c = P/[(n- 1)P] approaches the 
true age only in the limit Pq <C P. In the P-P diagram, inter- 
section of the spin-down trajectory P = Pn(P/Po) 2 " with an as- 
sumed spin-up line provides a recycled pulsar's initial period 
which, through equation (Q) yields the epoch since the end of 
spin-up. We apply this technique to estimate the ages of the 
known NS-NS systems; should accretion in a given system have 
halted before the equilibrium spin period was reached, or if the 
accretion rate was sub-Eddington, such an estimate becomes a 
strong upper bound on the true age. For the fiducial spin-up line 
and n = 3, spin-down trajectories for the known NS-NS systems 
are shown in Fig. 1, the P-P plane for recycled pulsars. 

The validity of post-accretion spin-down ages may be ques- 
tioned if there is no single spin-up line as given by equation ([[]). 
Evidence that accretion consistent with equation (jlj) is respon- 
sible for B1913+16-like systems (see Section |]) is provided by 
examination of the P-P diagram: all three of the known systems 
lie near the spin-up line. The pulsars in these binaries must have 
originated above and to the left of their present positions in this 
diagram, and in any case age estimates become increasingly in- 
sensitive to large values of the coefficient a (cf. Fig. 2). More 
importantly, a cannot be made arbitrarily big, as there are no 
known recycled pulsars above the spin-up line, excepting two 
pulsars in globular clusters (Section 5.1). We therefore adopt 



the value a = 1.6 x 10~ 15 s" 4 / 3 as a reasonable maximum to 
constrain ages at the upper end. 

We compute spin-down ages for several values of the brak- 
ing index, but for n ^ 3 the slope of the spin-up line differs from 
P 4 / 3 in a way that is difficult to quantify without a model for 
the braking torque produced by the magnetic field. We have ap- 
plied the model of Melatos (1997), which successfully explains 
the measured braking indices n < 3 of three young pulsars as 
a result of their magnetospheric structure, to objects near the 
spin-up line, and find that they should have n = 3 - e, where 
e is initially very small and decreases further as objects spin 
down. This result supports age estimates derived for spin-up 
lines oc P 4 / 3 and spin-down with n = 3. If magnetic fields of 
recycled pulsars evolve on timescales comparable to P/P, how- 
ever, n =/3 is possible (for specific models, see e.g. Ruderman, 
Zhu & Chen 1998); in such models, the spin-down time is al- 
tered, but not necessarily the equilibrium spin-up relation. 

2.3. Kinematic Ages 

A radio pulsar's distance z from the galactic plane together 
with its velocity derived from proper motion measurements 
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can be used to estimate its age through integration of the pul- 
sar's trajectory in the galactic gravitational potential. Sources 
of uncertainty in addition to measurement errors in distance 
and proper motion are the unknown line-of-sight velocities and 
non-zero z-heights at "birth" (the end of the spin-up phase 
for recycled pulsars). We calculate gravitational restoring 
forces towards the disk at the present-day galactocentric radii 
of B1534+12 and B1913+16 using the potential of Paczyriski 
(1990). 

3. RESULTS FOR THE KNOWN NS-NS BINARIES 

3.0.1. PSRB1913+16 

Useful kinematic age constraints can be derived for 
B1913+16 because it is at low galactic latitude and its distance 
and proper motion are well determined. Kinematic and spin- 
down age estimates are compared in Fig. 2. Trajectories in the 
top panel are constrained to cross the system's current z-height 
moving upward at the correct present-day velocity v z , they im- 
ply that the system is either leaving the galactic plane for the 
first time (r < 5 Myr), or has completed at least one-half cycle 
of its oscillation in z. In the latter case, i.e., for ages greater 
than a few tens of Myr, the system is likely to have sampled 
the galactic potential over a large range of galactocentric radii 
r g , so that the assumption of a constant r g breaks down — this 
will stretch or compress the plotted trajectories in a way that 
depends on the actual path taken by the system through the 
Galaxy, which is unknown. The unknown line-of-sight velocity 
v r also modifies these trajectories, but at the system's low galac- 
tic latitude, b = 2°, only a very large v r would be significant: to 
change v, by 10%, approximately the difference between tracks 
(a) and (b) in Fig. 2, requires v r ~ 400 km-s" 1 . 

The characteristic age of B 19 13+16, assuming a P corrected 
for galactic acceleration (Damour & Taylor 1991), is t c = 109 
Myr, while the combined kinematic and spin-evolution con- 
straints yield two possible solutions for the age: a young so- 
lution, t < 5 Myr, which is largely insensitive to uncertainties 
in the assumed distance, present-day v z , and braking index, and 
an old solution, 60-80 Myr, where the upper bound assumes 
n = 3 and q^1.6x 10~ 15 s" 4 / 3 . It is noteworthy that the old so- 
lutions provide an upper bound on the braking index, n < 3.5, 
unless a is very large. 

3.0.2. PSRB1 534+12 

The characteristic age of PSR B 1534+12 is 252 Myr, as- 
suming the acceleration correction for P given by Stairs et al. 
(1998). The pulsar's distance from the spin-up line (Fig. 3) im- 
plies t ^ 210 Myr, for n = 3 and a £ 1 .6 x 10" 15 s" 4 / 3 . The two 
estimates differ little because the condition Pq <C P is nearly 
met if the pulsar was born on the spin-up line. 

Based only on proper motion and distance, PSR B 1534+12 
appears to be moving toward the Galactic plane with transverse 
velocity 50 km-s" 1 from its current position at z — 800 pc. Be- 
cause the system is at Galactic latitude b = 48°, a modest radial 
velocity could reverse the sign of this apparent motion. Without 
more information about the present-day v, the system's kine- 
matic history is largely unconstrained. We have nevertheless 
derived the relevant motion in the Galactic potential under var- 
ious assumptions. In general, the minimum v z required at birth 
for the system to just reach its present z-height from a starting 
point zo = ± 200 pc is v z ~ 50 km-s . If we constrain the age 
to be roughly 200 Myr, allowed solutions span a range of initial 
velocities and maximum z-heights depending on the number of 



oscillation periods that have elapsed since the system's birth, 
but likely solutions require a birth v z between 100 and 200 km- 
s" 1 and Zmax of 2.0^1.5 kpc. Notably, the solution in which the 
pulsar has completed just less than half of its oscillation cycle 
and is returning to the plane for the first time requires a birth 
velocity v z ~ 275 km-s" 1 and z m ax = 13 kpc, which seem un- 
likely. If the age is in fact as large as we believe, the system has 
probably crossed the plane more than once. On the other hand, 
the kinematic age can be made arbitrarily small if the radial ve- 
locity is assumed to be large and the system is leaving the plane 
for the first time: an age ~ 10 Myr requires a v z at birth of 90 
km-s" 1 and a present-day v r ~ 140 km-s" 1 . 



3.0.3. PSRB2127+11C 

The characteristic age of B2 127+1 1C is t c = 97 Myr; the ap- 
parent P of this pulsar is likely contaminated by the system's 
unknown acceleration in the gravitational potential of M15 (but 
Phinney 1992 estimates the contamination to be at most 2%). 
Although the evolutionary history of B2 127+1 1C is murky, the 
pulsar's spin-down track in P and P probably originated some- 
where near the conventional spin-up line. We estimate the time 
since spin-up for this pulsar to be ~ 60 Myr (Fig. 3). 



PSR B 1 9 1 3+ 1 6 




Age (Myr) 

FIG. 2. — Kinematic and spin evolution constraints on the age of the 
PSR B 1913+16 system. The top panel shows trajectories perpendicu- 
lar to the galactic disk allowed by current measurement uncertainties 
in distance, d = 8.3 ±1.4 kpc, and proper motion (Damour & Taylor 
1991). The heavy horizontal line represents the system's present-day 
height z = 307 ± 52 pc above the disk. For each assumed distance, (a) 
6.9 kpc, (b) 8.3 kpc (emphasized), and (c) 9.7 kpc, four trajectories are 
plotted, one each for a birth z-height of ±200 pc and a present-day v z at 
either end of the la range, v z = +{133 ± 14)(rf/8.3 kpc) km-s^ 1 . Dashed 
curves are mirror-images through z = of the solid curves and show 
trajectories for negative-going birth z velocities; for ages between 50 
and 120 Myr, these are the relevant trajectories, as v z > today. In 
the lower panel, the dependence of spin-down age on the position of 
the effective spin-up "line" is shown for a variety of braking indices 
(curves labeled with n at top); the dotted line indicates the fiducial 
value of q used in our age estimates. The shaded region is kinemati- 
cally disallowed. 
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Time Since Leaving Spinup Line (Myr) 
Fig. 3. — Spin-down age constraints for PSRs B1534+12 and 
B2127+11C. See caption of Fig. 2 for details. 

4. THE FORMATION AND MERGER RATE OF NS-NS BINARIES 

4. 1 . Maximum Likelihood Birthrate Estimates 

For a given pulsar in a NS-NS binary, the birthrate of simi- 
lar objects may be estimated, following Phinney (1991), using 
AN = Vgai/VrfT, where V ga i is the total galactic volume in which 
NS-NS binaries are found, Vd is the luminosity- and density- 
weighted volume that pulsar surveys have sampled for objects 
with identical orbital period, and t is an appropriate source 
lifetime. The total birthrate follows by summing over objects, 
N = J2jANj. 

A more rigorous treatment, based on a likelihood function, 
uses all information in pulsar surveys, including non-detections 
in much of the survey volume as well as detected pulsars. Al- 
though the scarcity of known binary pulsars would frustrate any 
attempt to constrain meaningfully the many important param- 
eters that characterize the Galactic population of NS-NS bi- 
naries, we can nevertheless use the formalism to understand 
the proper interpretation of the parameters that enter into any 
birthrate estimate. We find that a maximum likelihood estimate 
for the binary neutron star birthrate N is 



N-- 



N d 



(fr] 



(5) 



exp 



where the number of detected binary pulsars is Nj and {fT) exp 
is the joint expectation value of the fraction of all existing bi- 
nary pulsars detected times source lifetime, summed over sur- 
veys. The denominator involves an integral of the number den- 
sity of binary pulsar systems (normalized to one) times their av- 
erage lifetime convolved with observational selection functions 
that limit the spatial volume and range of lifetimes amenable 
to detection. At a more fundamental level, the necessary inte- 
gral involves physical distributions of the positions, center of 
mass speeds, orbital elements and luminosities of binary pul- 
sars at birth, and also requires physical models for the evolution 
of these distributions along with detailed understanding of the 
phenomenological limitations of searches for such systems. 



Ideally, the physical distributions involved in computing 
(f T )exp should be determined from the available data, but, 
since very few NS-NS binaries are known, we, like others be- 
fore us, are forced to appeal to prior knowledge or theoreti- 
cal expectations about the binary populations instead. More- 
over, we approximate the expectation value for each survey by 
(fT)exp ~ {f)ex P {T)exp, where (f) ap is the fraction of all NS- 
NS binaries we expect to be observable, and (r) eX p is the de- 
tectable lifetime of such systems (see, e.g., Van den Heuvel & 
Lorimer). We can also approximate 



{f)exp 



fdVd 



(6) 



where fa < 1 is a measure of how representative the detected 
binary pulsars are of the overall Galactic population. If the NS- 
NS binaries found so far are typical of the overall population 
then fa « 1; if not, then fa could be small (e.g., a substan- 
tial fraction of binary pulsars could have orbital periods short 
enough that detection is unlikely; Johnston & Kulkarni 1991, 
Bloom, Sigurdsson & Pols 1998). For fa = 1, we find a birthrate 
estimate from a single detected binary pulsar similar to Phin- 
ney's (1991) AN, except that it is now clear that (i) there is 
an underlying assumption that the detected system is typical of 
the population and (ii) the appropriate lifetime to use is the ex- 
pected time over which the pulsar should be detectable. Finally, 
in the absence of detailed knowledge of the spatial distribution 
of NS-NS binaries, we shall approximate (following the lan- 
guage of Cordes & Chernoff 1997) the "detection volume" Vd 
by the "sampled volume" V s , which is the sensitivity-limited 
region searched by pulsar surveys. This approximation is con- 
servative in that V s overestimates Vd (see below). 

Apart from fa, a major uncertainty in any birthrate estimate 
is the value of (T) exp . Below, we estimate this parameter by 
computing the expected detectable lifetimes for B 19 13+ 16 and 
B1534+12 separately. How to translate these values into an es- 
timate of (r)«xp is somewhat subjective, as it depends on one's 
prior notions of which, if either, of them is "more typical." This 
subjectivity pervades previous estimates of N as well as our 
own: all are proportional to nonrigorous estimates of (t)~ xp . As 

should be evident from equation (|j), evaluating N as the sum of 
birthrate estimates for individual binary pulsar detections rep- 
resents a particular choice for (r)~ xp . 

In a manner similar to Phinney (1991), the birth and merger 
rates of NS-NS systems have been estimated by Narayan, Pi- 
ran & Shemi (1991) and Curran & Lorimer (1995) under, ef- 
fectively, different assumptions for V ga i and V s . In each case, 
the lifetimes of individual systems were taken to be the sum 
of the pulsar's characteristic age and the binary's coalescence 
time. Curran & Lorimer's (1995; hereafter CL95) estimate 
has twice been revised. Van den Heuvel & Lorimer (1996; 
hereafter VL96) argued that the effective "detectable lifetime" 
(rather than the full lifetime culminating in merger) should be 
used in birthrate calculations; they estimated radio luminos- 
ity decay timescales for each pulsar and proposed a merger 
rate nearly three times that of CL95. Stairs et al. (1998; here- 
after SAC+98) determined that the distance to PSR B 1534+12 
is larger than previously believed, implying that V s had been 
underestimated by CL95. Their correction reduced Van den 
Heuvel & Lorimer's (1996) estimate by a factor of 2.5, bring- 
ing it effectively back to CL95's value. These revisions are 
summarized in Table 1 . 



5 



Table 1 

Recent (post-1995) birthrate estimates for NS-NS binaries. For each entry, the first item in the pair refers to 

B1534+12, the second to B1913+16. 





Scale factor, £ 


Age, r 


Remaining detectable 


Detectable lifetime, 


Individual rate, 






(Myr) 


epoch, t> (Myr) 


T = r + T r (Myr) 


zvv = £/r(yr I ) 


CL95 


160, 20 


250, 109 


2700, 300 


2950, 409 


5.42 x 10" 8 ,4.89x 10" 8 


VL96 


160, 20 


250, 109 


510, 220 


760, 329 


2.10 x 10" 7 , 6.08 x 10" 8 


SAC+98 


37.8,20 


250, 109 


510, 220 


760, 329 


4.97 x 10" 8 , 6.08 x 10" 8 


our revision 


37.8,20 


200, 65 


510, 180 


710, 245 


5.32 x 10" 8 , 8.16 xlO" 8 


present work" 


V ga i/V s = 107.1 


200, 65 


510, 180 


r = 500 


V g ai/V s r = 2.14x 10"' 



"This entry summarizes our estimate N^$, equation (^|), which is based on different assumptions for computing galactic scale factors 
than the first four entries (see text). 



z distribution of old NS which, for isolated pulsars, includes a 
large fraction that oscillate to z max ~ 10 kpc. We argue that the 
space velocities of NS-NS binaries are sufficenfly large (up to 
~ 200 km-s" 1 ) that both the radial and z distributions should be 
extended, to roughly a 10 kpc characteristic radius and h z m 5 
kpc. This yields V ga i w 3000 kpc 3 . 

Estimates of V s are highly sensitive to modeling of selection 
effects in surveys that have or could have found the NS-NS bi- 
naries. The volume sampled (per unit solid angle) in a survey 
is (e.g., Cordes & Chernoff 1997) 



dVs 



= / dL p fi p (Lp) 



>V2 



dDD 



(7) 



We first report here our own revision consistent with the 
methods and assumptions of CL95. We then evaluate the as- 
sumptions entering into previous V ea i and V s calculations and 
report a birthrate reflecting recently completed radio pulsar sur- 
veys. 

4.2. Effects of Revised Lifetimes on Birthrate Estimates 

As argued by Van den Heuvel & Lorimer (1996), and as 
emphasized above, we need the detectable lifetime T of the 
sources to estimate their birthrate. VL96 considered radio lu- 
minosity decay and showed that both B 1534+12 and B 19 13+ 16 
will likely become undetectable before they merge with their 
companions. However, detection is biased against finding recy- 
cled pulsars in short-period binaries; in present surveys, neither 
pulsar might have been found if their binary periods were < 0.5 
times their present values. The gravitational radiation lifetimes 
of B1534+12 and B1913+16 are 2.7 x 10 9 yrs and 3.0 x 10 8 
yrs, respectively. For B 1534+12, the orbital period will halve 
in about 2.2 x 10 9 yrs and for B1913+16, in about 1.8 x 10 8 
yrs. The latter result implies that the detectable lifetime of 
B 19 13+ 16 is limited by orbital selection rather than by lumi- 
nosity decay. As shown in Table 1, then, our revision to CL95's 
estimate consists of adopting (1) the Stairs et al. (1998) distance 
to B1534+12, (2) our preferred age estimates for B1913+16, 
65 Myr, and B 1534+12, 200 Myr, and (3) retaining Van den 
Heuvel & Lorimer's (1996) luminosity decay for B 1534+12, 
but assuming B 19 13+16 will remain detectable for 180 Myr as 
a result of observational selection due to the rapidly decaying 
orbit 1 . We thus derive a birthrate, uncorrected for beaming or a 
luminosity function to account for dim objects, of 1.35 x 10~ 7 
yr" 1 by summing the individual rates of the two objects, a 23% 
increase over the latest revision. We note that CL95 suggest a 
factor of 10 increase to account for low-luminosity objects, and 
a further factor of three increase for beaming. 

4.3. Effects of Revised Galactic Volumes on Birthrate 
Estimates 

Values of V ga i used previously range from ~ 100-900 kpc 3 ; 
the lowest V ga i (CL95) attributes a Population I radial distri- 
bution and a scale height h z = 1 kpc to NS-NS binaries. The 
kinematic "ages" of binary pulsars are directly related to their 
scale heights. Previously, a range h z = 0.5-1 kpc has been used 
to estimate binary pulsar birthrates, but this does not typify the 

'While relativistic precession of the pulsar spin axis may render B 19 13+16 undetectable on much shorter timescales (Kramer 1998, and references therein), for a 
population of objects only the static beaming properties represented here by the beaming fraction f}, determine the typical detectable lifetime. 



where Z) max = (L p /S m i n ) is the maximum distance sampled, 
S m in is the minimum detectable flux density, and fi_ is the lu- 
minosity function (with unit area) of the "pseudo luminosity" 
L p = D 2 S4oo, where S400 is the 400 MHz flux density tabulated 
in Taylor, Manchester & Lyne (1993). Equation (^) provides a 
good estimate of the detection volume Vd when the scale height 
h z of the NS-NS population is high enough, that is, when even 
the brightest pulsars in the population cannot be detected be- 
yond a vertical distance h z above the location of the solar sys- 
tem in the Galactic disk; for smaller h z , equation (Q) overesti- 
mates Vd- Previous work indicates fi p oc L~ 2 over several orders 
of magnitude of L p , for both isolated pulsars (Lyne et al. 1998) 
and millisecond pulsars (Cordes & Chernoff 1997). The lower 
cutoff is taken to be L pi ~ 1 mJy kpc 2 . For the upper cutoff, we 
use the value for B1913+16, L ;; , ~ 200 mJykpc 2 , as its pulse 
shape suggests that the radio beam points directly at us, yield- 
ing a large, perhaps maximal luminosity for its P and P. This 
luminosity function does not account for pulsars beamed away 
from us completely, so our rates still must be adjusted upward 
for beaming. 

The minimum detectable flux density, S m \ n , is generally a 
strong function of P, DM, P or b, and direction, as well as of in- 
strumental parameters. For long spin- and orbital-period pul- 
sars, S m i n asymptotes to a constant for low DM. Using the 
asymptotic 5 m ; n ~ 0.8 mJy representative of recent surveys at 
Arecibo, we obtain 



dV s 
dTL 



< -<T 3/2 L 

<■ 3 ■V- L 



Pi L P2 



1/2 . 



10- 2 - 4 kpc 3 deg- 2 ^ 2 . 8 L Pl , 1 L 



1/2 

P2.200' 

(8) 
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where S min = 0.85^,0.8 mJy, L PI = IL PiA mJykpc 2 , and L P2 = 
200L P2 ,2oo mJy kpc 2 . Arecibo surveys have covered Q s « 5000 
deg 2 (Camilo 1998) in recent years; surveys at Parkes and Jo- 
drell have covered more solid angle but with S m i n ~ 5 mJy. To- 
gether, these imply V s <, 28 kpc 3 . We stress that this value may 
be a significant overestimate, for two reasons. First, the asymp- 
totic flux S m \ n is not appropriate for rapidly rotating, highly dis- 
persed pulsars, especially at low galactic latitudes. Second, the 
true low-end cutoff of the radio luminosity distribution could 
occur at luminosities much lower than 1 mJykpc 2 — Lyne et al. 
(1998) show that the luminosity function for non-recycled pul- 
sars begins to turn over for L p < 10 mJykpc 2 , but that there is 

no evidence for a turnover at 1 mJykpc 2 for spun-up objects. 

Using our estimates of V ga i and V s , and expressing (T) exp in 
units of the mean of the detectable lifetimes of B 19 13+16 and 
B 1534+12, we obtain a rate 

NS y boOOkpcVV {r)«p J\ V s ) Jh 

(9) 

where fi, is the beaming fraction, under the assumption that 
these binaries are typical of the Galactic population of NS-NS 
systems. This result is nearly a factor of 7 smaller than our re- 
vision of the CL95 estimate corrected for their suggested factor 
of 10 contribution from low-luminosity objects, lO" 5,9 /^ 1 yr" 1 . 
We attribute this difference to the increased galactic volume re- 
cently surveyed for rapidly rotating pulsars as represented here 
by V s , and our inclusion of the radio luminosity function in esti- 
mating V s instead of the somewhat arbitrary multiplying factor 
ofl0ofCL95. 

Aside from the beaming factor fy, equation (^) is a lower 
bound to the actual birthrate for three reasons. First, the 
birthrate function for binary neutron stars could be skewed to- 
ward short orbital periods (Bloom, Sigurdsson & Pols 1998) 
where detection becomes unlikely (Johnston & Kulkarni 1991). 
If this is true, then fj is small, and the birthrate is larger than 
iv£,s by a factor of f~ t l (see eqs. [^] and [^]). Second, V s over- 
estimates Vd'. for high-sensitivity surveys such as those carried 
out at Arecibo, sky coverage at high latitudes samples a large 
volume at z-heights beyond the likely scale height of NS-NS 
binaries. Third, the parameters entering into our estimate of V s , 
equation (|j), are sufficiently uncertain that the sampled volume 
itself is likely overestimated. We have probably overestimated 
(r) exp as well. Taken together, these considerations suggest that 
nLI may underestimate the true rate by an order of magnitude 
or more. 

Previous estimates of the binary pulsar birthrate, includ- 
ing our N^l, presume that their lifetimes are limited by ei- 
ther spiral-in or spin-down. Ideally, one should dispose with 
this prior assumption and determine (r) eAp from the data, if 
possible. We note that kinematical considerations imply that 
B 19 13+16 may be ~ 10 Myr old, and a similarly young age 
cannot be ruled out for B1534+12. If a significant fraction 
of recycled pulsars discovered in future surveys prove to be 
young (age ~ 10 Myr), but in binary systems with long gravi- 
tational radiation and spindown lifetimes (~ 100-1000 Myr), 
then we would be forced to reduce our estimate of the de- 
tectable lifetime r. (This would also imply that something other 
than gravitational radiation or spindown regulates the timescale 
over which binary pulsars are observable.) Such a short lifetime 

2 The NS-NS system B2303+46, which is known through the non-recycled r 



would require a large birthrate, NBl ~ 10~ 5 l) f^ 1 yr -1 , for such 
systems, which could dominate the rate of neutron star mergers 
if they are common. As long as there is not a statistically sig- 
nificant population of anomalously young binary pulsars, this 
large birthrate is not required. Given that only a few pulsars are 
known to exist in tight binaries, however, we cannot rule out the 
existence of short-lived binary pulsars that are born frequently. 

Both of our birthrate estimates and N^l can easily take 
on values larger than 10~ 5 yr -1 when provision is made for 
beaming and the large uncertainties in V ga i, V s , and Vd. This pos- 
sibility prompts us to re-examine the upper bound of 10~ 4,9 yr" 1 
proposed by Bailes (1996) based on the absence of any known 
"normal" (i.e., non-recycled, long-period, high magnetic-field) 
pulsars in a close NS-NS binary. We begin by defining the fol- 
lowing quantities: 

• N, the birthrate of NS-NS binaries that will merge in a 
Hubble time, 

• T„b, the detectable lifetime of the non-recycled neutron 
star in such a binary, 

• f V b, the fraction surveyed to date of the galactic volume 
occupied by such binaries (= Vd/V ga i), 

• /orb, the fraction of such binaries that are detectable de- 
spite orbital smearing of the non-recycled pulsar's pulsa- 
tions, 

• N n , the birthrate of (normal) isolated pulsars, 

• T„, the detectable lifetime of a normal pulsar, 

• /„„, the fraction surveyed to date of the galactic volume 
occcupied by normal pulsars, 

• N„, the number of non-recycled pulsars known to date. 

As no normal pulsars have been found in NS-NS systems ex- 
pected to merge in a Hubble time, 2 we have 

NxT nb x f vb x / orb x f b < 1. (10) 

Similarly, for the normal, isolated pulsar population, 

N n xT n x f m x f b =N„. (11) 

Taking the ratio, and assuming T„b = T n , we get 

N<^(-^-)Nn. (12) 

N n V/orb/W 

The quantity in parentheses, which was not considered by 
Bailes (1996) explicitly, can be less than or greater than one, 
but is probably close to unity. We believe the scale height of 
binaries is somewhat smaller than that of the isolated normal 
pulsars, so f v b ;> f vn is likely (i.e., a larger fraction of the galac- 
tic volume of binaries has been searched). Because of orbital 
acceleration, / or b is less than one, but represents a less severe 
selection effect for slow pulsars than for recycled objects. 

Bailes (1996) assumed N„ = 650, and adopted the largest al- 
lowed birthrate, B„ = 1/125 yr" 1 , from the analysis by Lorimer 
et al. (1993) of the population of normal, isolated pulsars. A 
similar recent effort by Lyne et al. (1998) based on a larger sam- 
ple of objects finds a maximum birthrate of one every 60 yrs. 

star, has P or (, = 12.3 days and will not coalesce in less than a Hubble time. 
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Furthermore, the number of cataloged normal pulsars with mea- 
sured spin-down rates, i.e., for which non-membership in a NS- 
NS binary has been firmly established, is 501. Assuming unity 
for the factor in parentheses in equation (|T2|), an upper bound 
on the birthrate of close NS-NS binaries is then N < 10~ 4 5 yr -1 , 
independent of beaming. This rate should be corrected upward 
by a statistical "Poisson" factor and by the factor / or b; we there- 
fore believe that a robust upper limit is 

A r <10 _4 yr- 1 . (13) 

We point out, however, that the rate N$,£ need not be subject 
to this bound, as the assumption T„b = T„ might not hold if the 
detectable lifetimes of NS-NS binaries were somehow not gov- 
erned by the spin-down or gravitational radiation timescales. 

Finally, we refer, for completeness, to the recent work of 
Rosswog et al. (1998) who find that a NS-NS merger rate of 
10~ 5 5 -10~ 4 5 yr" 1 could acount for the apparent deficit of r- 
process nuclei. 

5. ACCRETION MECHANISMS AND THE DISTRIBUTION OF 
RECYCLED PULSARS IN P AND P 

Surveys over the last 16 years have been sensitive to objects 
with P J> 1 ms to varying degrees. Certain portions of the P-P 
diagram have become populated as a result, while others have 
remained empty. We argue that three of these unpopulated re- 
gions, labelled A, B, and C in Fig. 1, provide insights into the 
evolution of binary systems containing neutron stars. In partic- 
ular, there is a dearth of pulsars in region B where we would 
expect to find inspiralling NS-NS binaries. We assess the sig- 
nificance of this and other deficits below. 

5.1. Region A 

There are no known selection effects to prevent detection of 
pulsars in region A. objects with spin-down rates higher than 
the low-field millisecond pulsars would likely be brighter and 
so more easily detectable. Still, no recycled pulsars in the galac- 
tic disk are known above the spin-up line. We interpret this as a 
qualitative validation of the Ghosh-Lamb spin-up line. We dis- 
count pulsars in globular clusters, here and in the discussions to 
follow, because of their uncertain spin-down rates and compli- 
cated evolutionary histories. 

5.2. Region B 

Region B is below and to the right of the known NS-NS sys- 
tems, for 50ms <; P <^ 200ms (Fig. 1). A deficit of known pul- 
sars here can occur for three reasons, all of which likely play 
some role: 

1. NS-NS binaries born near the spin-up line will coalesce 
before reaching region B. In Fig. 1, the spin-down tracks 
plotted for the three NS-NS systems end at lower right 
at the expected coalescence time for each system. Only 
B 1534+12 penetrates into region B before merging. 

2. For systems like B 1534+12, inspiral results in short or- 
bital periods, and high orbital velocities, within region B, 
and such pulsars have been selected against. 

3. As the pulsars in these systems age, their radio luminosi- 
ties probably decay, making them more difficult to detect 
(Van den Heuvel & Lorimer 1996). 

3 We predict that should the NS companion to B2303+46, a 1 s pulsar in a 12 
alongside the similar NS-NS system J 15 18+4904. 



Evidently, any systems in B are selected against in pulsar sur- 
veys, or few live long enough to reach or spend sufficient time 
in the region. Apparently, no new systems like B 19 13+16 are 
being formed far from the spin-up line. The other known pul- 
sars in the vicinity of region B are consistent with this hypoth- 
esis. The 0.8 Mq white-dwarf companion to PSR B0655+64, 
the "failed NS-NS" system, has cooling age ~ 4 Gyr (Hansen 
& Phinney 1998), which suggests the pulsar emerged from the 
accretion phase near the spin-up line with Pq ~ 20 ms. Also, 
J2235+1506 and two objects with P > 200 ms are solitary and 
so could not have escaped detection through merger or orbital 
acceleration. The suggestion of Camilo, Nice & Taylor (1993) 
that J2235+1506 is the result of a disrupted binary is consistent 
with this view. 

The deficit of pulsars in region B supports our use of the 
spin-up line as a zero-age starting point, and is consistent with 
the existence of a population of NS-NS systems with periods 
shorter than 8 hours which has remained undetected because 
of sensitivity limits and observational selection against binary 
pulsars with short orbital periods and massive companions. 

5.3. Region C 

Pulsars should already be known in region C, periods 5- 
30 ms and P > 10~ 19 s-s" 1 , where they would have spin- 
down luminosities 1 .5-2 orders of magnitude higher than, say, 
J0621+1002. We assert that the evolutionary path for producing 
radio pulsars in region C is disallowed or rarely followed. 

A clue to the absence of pulsars in C may lie in the ap- 
parent bimodality in the spin parameters of the massive bi- 
naries (Fig. 1). The "B1913" family (B1913+16, B1534+12, 
B2 127+1 1C, and B0655+64) consists of pulsars in short-period 
orbits with NS (or white-dwarf) companions; the pulsars had 
P{) > 20 ms at the end of the spin-up phase. The "J0621" fam- 
ily 00621+1002, J1518+4904, J2145-0750, and J1022+1001) 
consists of pulsars in ~ 8-day orbits with massive white-dwarf 
(or NS) companions. The division between the two families is 
strikingly clean — the systems with discrepant companion types 
(B0655+64 and J15 18+4904) nonetheless have orbital periods 
consistent with their respective family members. 3 

The shorter average spin period of the J0621 family implies 
that spin-up in the wide-orbit systems must have been at least 
as effective as in the close NS-NS binaries, but produced fi- 
nal magnetic fields smaller by more than an order of magni- 
tude. These facts suggest that different accretion and spin-up 
processes were responsible: we believe that the B1913 fam- 
ily was spun up through disk accretion and the J0621 family 
predominantly during a wind-driven accretion phase (e.g., Van 
den Heuvel 1994 presents a spin-up scenario for J2145-0750 
involving wind accretion). It is generally believed that disk ac- 
cretion drives monotonic spin-up (until the spin period ~ P eq ) 
once the donor star overflows its Roche lobe, whereas wind ac- 
cretion produces a non-monotonic random walk in the angular 
momentum imparted to the NS. To achieve the same final spin 
period, then, much more mass must be transferred, at lower ac- 
cretion rates, in wind accretion, resulting in smaller magnetic 
fields for the wide-orbit binaries; this is consistent with viable 
models of field attenuation during accretion (Konar & Bhat- 
tacharya 1998). (Wijers 1998 refutes the hypothesis that field 
loss scales with accreted mass, but without considering non- 
monotonic spinup.) Although recent observations show that 

eccentric orbit, ever be detected as a pulsar, it will lie among the J062 1 family, 
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disk accretion is also not entirely monotonic (Bildsten et al. 
1997; Nelson et al. 1997), and certain numerical simulations 
of wind accretion question the existence of angular momentum 
"flip-flops" (e.g., Font & Ibanez 1998), frequent or deep torque 
reversals during spin-up in wide binaries can account for the 
two distinct families of massive binary pulsars without the need 
for abandoning field decay that scales with accreted mass. 

Given the proximity of region C to the spin-up line and the 
implied high magnetic fields relative to J062 1+1002, our argu- 
ments suggest that the missing NS undergo accretion in close 
systems. We speculate, then, that the initial conditions (orbital 
separation and companion type) needed to produce pulsars in 
region C lead instead to accretion ending with collapse of the 
NS to a black hole or merger with its companion in a common- 
envelope phase. In detail, the outcome would depend on the 
maximum mass of neutron stars and the efficiency of the accre- 
tion (see, e.g., Brown 1995 and Fryer, Benz, & Herant 1996). 

Detecting pulsars in region C should also be made difficult by 
the same pulse-smearing effects that would prevent detection 
of pulsars like B 19 13+ 16 in binary orbits with periods smaller 
than about four hours (Johnston & Kulkarni 1991; Section 4 
above). The minimum detectable spin period for a given or- 
bital period is P m { a oc P ^ 3 , so it is unlikely that pulsars spin- 
ning faster than P ~ 20 ms could be detected in binaries with 
f 01 b ~ 8 hrs. Conversely, pulsars with P ~ 20 ms could only 
be found in binaries with orbital periods longer than about 8 
hours. This means that recycled pulsars in binaries like those 
that contain B 1534+12, B 19 13+ 16, and B2 127+1 1C could re- 
main undiscovered throughout region C. Moreover, the absence 
of pulsars in somewhat longer period binaries (e.g. f or b ~ 0.5-1 
day) in region C suggests that there is a fairly sharp upper cutoff 
to the distribution of periods of close binary pulsar systems, the 
extended B1913 family. This supports the idea that the much 
longer period binaries in the J0621 family have a completely 
different physical origin. 

5.4. Discussion 

The explanation commonly put forth for millisecond pulsars 
with t c exceeding a Hubble time is that these objects were born 
with spin periods close to their present values: spin-up end- 
points are assumed to adhere to equation ([!]), but with very 
low accretion rates, m ~ 10~ 2 (e.g., Hansen & Phinney 1998). 
We suggest two additional ways in which spin-up endpoints far 
from the canonical spin-up line may be produced: (1) for short- 
period binaries (P or b < 2 d), pulsars emerge on one of the P~ s 
spin-up lines if accretion ceases before the equilibrium period 
is reached; and (2) for long-period binaries, wind accretion re- 
sults in low-field radio pulsars with as yet undetermined initial 
periods. The ages of such systems are therefore uncertain. It 
would seem that near-Eddington disk accretion given sufficient 
time to reach equilibrium rotation and ending in a pulsar is rare, 
as is indicated by the paucity of radio pulsars near the spin-up 



line for P < 30 ms. In addition, if prolonged mass transfer with 
minimal spin-up during the wind accretion phase is an impor- 
tant part of the evolution of wide-orbit millisecond pulsars, as 
we suggest, the masses of these objects may be considerably 
higher than previously believed. Future high-precision mass 
measurements for such systems may therefore provide impor- 
tant constraints on the maximum allowed neutron star mass. 

6. CONCLUSION 

By estimating the detectable lifetimes of PSRs B1534+12 
and B 19 13+ 16, we have revised previous estimates of the galac- 
tic birthrate of close double neutron star binaries. In steady 
state, this would also be the coalescence rate for these relatively 
short-lived binaries. We find, adopting the galactic "scale fac- 
tors" of Curran & Lorimer (1995), that the birthrate of such 
systems, uncorrected for a radio pulsar luminosity function and 
beaming, is 10~ 6 9 yr" 1 . Our own estimate of the birthrate, inde- 
pendent of Curran & Lorimer's (1995) methodology, is based 
on a maximum likelihood analysis that clarifies the assumptions 
entering into all birthrate estimates, including our own. We 
find, given available knowledge of the ages, systemic velocities, 
and orbital properties of the two known close NS-NS systems in 
the galactic disk, along with an estimate of the galactic volume 
surveyed to date, that the birthrate (and hence merger rate, for a 
steady-state distribution) is J£ lO -6 ' 7 /^" 1 /^ 1 y r_1 > where fa is the 
beaming factor, and fa < 1 represents the degree to which the 
two known objects are typical of the entire population of NS- 
NS binaries. Important parameters entering into this estimate, 
e.g., the lower cutoff of the luminosity function for recycled 
pulsars, are sufficiently uncertain, however, that the true rate is 
probably higher by an order of magnitude or more; we therefore 
consider our result to be a realistic lower bound on the true rate. 
Furthermore, we find that kinematic and spin-down age con- 
straints for PSRs B1534+12 and B1913+16 allow both systems 
to be unexpectedly young — if we dispose of prior prejudices 
about the detectable lifetimes of these systems, their birthrates 
could be substantially higher, of order 10~ 5 fa l fa l yr" 1 . We 
have revised the upper bound on the NS-NS binary birthrate 
proposed by Bailes (1996) based on the absence of any known 
non-recycled pulsar in a NS-NS binary, finding that the po- 
tentially large birthrates allowed by our analysis are consistent 
with a revised upper bound at 10~ 4 yr" 1 . We therefore bracket 
the merger rate of NS-NS binaries to the range ^0~ 67 fa 1 fa 1 
yr-^JV^lO-V" 1 . 

Our considerations of V ga i for NS-NS binaries suggest that 
searches for such binaries out of the Galactic plane will be more 
productive than those in the plane; independent of how many 
are found, high-latitude searches will allow determination of 
the thickness of the population in z. 

This work was supported by the NSF and NASA grants 
AST95-28394 and NAG5-3097. 
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